Macroscopic Quantum Coherence in a repulsive Bose-Einstein condensate. 
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We consider a Bose-Einstein bicondensate (BEC) of ^'^ Rb, trapped in two different internal levels, 
in a situation where the density undergoes a symmetry breaking in momentum space. This occurs 
for a suitable number of condensed atoms within a double well dispersion curve, obtained by Raman 
coupling two internal states with two tilted and detuned light fields. Evidence of bistability results 
from the Gross-Pitaevskii equation. By second quantization, we evaluate the tunneling rate between 
the two asymmetric states; the effects of losses on coherence are also considered. 



The transition from a superfiuid to a Mott insulator 
has been recently demonstrated for a Bose-Einstein con- 
densate (BEC), made of atoms with repulsive mutual in- 
teractions in a lattice potential Such a phenomenon 
goes beyond the mean field approximation and its expla- 
nation requires to take into account the quantum fluc- 
tuations of the matter field. When the quantum tun- 
neling between adjacent sites dominates the interaction 
energy, the probability distribution for the atomic occu- 
pation of a single site is Poissonian. In the opposite case 
the minimum energy is obtained reducing the quantum 
fluctuation of the local occupation number. 

In a previous paper we have studied a similar problem, 
but with only two wells and an attractive interaction ; 
In this case the opposite effect occurs, that is, raising the 
atom number, i.e. the interaction energy, the minimum 
energy state is obtained increasing the atomic fluctua- 
tion in each site. This is indirectly demonstrated by the 
numerical observation of the symmetry breaking at a crit- 
ical number of atoms. Approaching the threshold value, 
the quantum fluctuations increase, whereas above that 
value they blow up and a new minimum quantum state 
appears. The associated probability distribution of the 
condensate barycenter displays two peaks, that can be 
considered as the dead and alive states of a Schrodinger 
cat (SC), whose coherent superposition is called Macro- 
scopic Quantum Coherence (MQC). This term was in- 
troduced to describe the coherent superposition of two 
macroscopically distinct quantum states that differ for 
the value of a collective variable Q . The phenomenon is 
observable only for attractive interactions, that tend to 
cluster the atoms in one of the two wells. By contrast, a 
repulsive interaction tends to reduce the quantum fluc- 
tuations and to distribute the same number of atoms in 
each well, as observed in Ref. |l|. MQC has been observed 
with trapped ions [| and microwave fields in high-Q cav- 

ity @. 

In this work, we discuss the feasibility of MQC in a 
BEC of mutually repulsive atoms. A repulsive interac- 
tion acts in the momentum space as an attractive inter- 
action, therefore we expect that for a double well dis- 
persion curve the symmetry breaking occurs in the re- 
ciprocal space. Such a dispersion curve can be obtained 
by two detuned and tilted light fields, that transfer a 



net momentum to atoms as they jump from an inter- 
nal state to another. We study the problem by finding 
the stationary solutions of two coupled Gross-Pitaevskii 
equations (CP) discretized over a space lattice. The sys- 
tem undergoes symmetry breaking in momentum space 
for a suitable number of atoms and two new stationary 
states are created. We then introduce a quantum two 
mode model, with the two modes chosen in such a way 
as to reproduce the stationary solutions of CP, and eval- 
uate the quantum fiuctuation-mediated tunneling rate 
between the two asymmetrical states. If the coupling 
with the environment is negligible, MQC occurs between 
these states. A Raman scheme for creating a superposi- 
tion state with two Rb condensates in different internal 
quantum levels has already been discussed in Refs. ^, |^; 
however both papers limit themselves to co-propagating 
light beams, and this implies applicability problems, as 
discussed in Ref. ^. 

Here we refer to atoms, but our numerical results 
apply also to '^^Na, if some parameters are rescaled. In 
a previous work Q we considered atoms in two different 
hyperfine levels {F — 1, mp = —1 and F = 2, mp = 1), 
however the associated depletion rate |^ is too high for 
our purposes. Here we consider condensate atoms that 
are optically trapped in the two Zeeman levels F = 1, 
mp = —1 and F = 1, mp = 1. An all optical con- 
densation has been reported in Ref. alternatively, the 
condensate can be created with a magnetic confinement 
and transferred into an optical trap . A homogeneous 
magnetic field has to be applied to remove the energy 
degeneracy. These levels are quasi-resonantly coupled by 
means of two Raman fields L and R. We call ipo and V'l 
the fields associated with the mp = ~l and mp — 1 lev- 
els, respectively. Furthermore we call ^'2 the upper state 
of the Di transition. 

The starting equations are 

(1) 

ihij;2 = huj2il^2 + [nEl{t)e'^^^-^-'^^'HQ + {L ^ R)^i] (2) 

where = -ffi+.goo|'/'oj^+.9oi|'0iP, Hi ^ i?;+£ii|7/'ip-t- 
5io|V'oP and gij = Aufi aij/m. Oij ^ 5.5nm |I3] are the 
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s-wave scattering lengths between atoms in i and j lev- 
els. We have called Hi = -(?i^/2m)V^ + V the one 
atom term of the Hamiltonian, where V is the trapping 
potential. The field amplitudes E^ ji are rescaled in or- 
der to be expressed in frequency units. They are thus 
the Rabi frequencies of the one-photon transition. huj2 
is the energy of the upper state of the one photon tran- 
sition; hujQ^i are the energies of the mp — —1, 1 levels, 
respectively. We set luq = 0. In the adiabatic approx- 
imation, '02 can be expressed in terms of ipi and ^ 
^2 = -[£;£e*(^^-^-"^*)Vo + £^Re'(^«-^-'^«*Vi]/A, where 
A = LU2 — lol- Thus we have two closed equations for 0o 
and 01. 

We introduce the gauge transformation 



V'o,: 



where kfi = kj^ — kji and = lol 
equations of motions become 



-^*±^*)V'o,i, (3) 
LOj^. As a result, the 



— 7y V'oa (4) 

2m 



Here, f2 = is the two photon Rabi frequency, taken 

for simplicity as time independent and real, and the fre- 
quency S is given hy S = uJi — uJd + (|£^Lp — |^^i?P)/A. 
We assume \El\'^ — |i?_R,|^ and uJd = t^i, hence 5 — 0. 

If the number of atoms is sufficiently small, we can 
neglect the nonlinear terms. Furthermore, let us ini- 
tially consider a spatially homogeneous condensate (no 
trap potential). As a consequence, Eqs. reduce to 
two linear equations with constant coefficients, and the 
eigenvalue problem in the reciprocal space is ruled by 
two linear algebraic equations for the transformed fields 
4>o{k) and 4>i{k). The momenta of the atoms in the two 
levels mp — —1 and mp — 1 are respectively h{k — kd/2) 
and h{k + kd/2). 

Solving the eigenvalue problem we find the two dis- 
persion curves hu>{k) = ti^k'^/(2m) ± {[h'^kdk/{2m) + 
M/2]2 + ?i2|rj|2}i/2^ plotted in Fig. §i for 17 = 200s-^ 

and \kd\ = 2tt ■ 4.5 • lO^m"^. As shown in the figure, 
the electromagnetic coupling modifies the parabolic dis- 
persion curves associated with the two hyperfine levels, 
lifting the degeneracy at their intersection point. The 
energy gap for S — is 2?i|f2|. By varying S one can rise 
or lower the energy separation between the two minima. 
For S — the two minima have the same energy. In- 
troducing the harmonic trap potential, the ground state 
has no longer a definite momentum. Furthermore, the 
two wells of the dispersion curve are equally populated 
by quantum tunneling (Fig. |^b, solid line). 

When the number of atoms is sufficiently high, a sym- 
metry breaking occurs because of the atomic interactions. 
We report in Fig. [l|3 the distributions |0oP and |0ip for 
(5 = and for three different values of the number of 
atoms N. 0o,i sue the Fourier transforms of the ground 
state solution (00, 0i) of Eqs. (Q) for a spherical trap po- 
tential V corresponding to equal longitudinal (wj_) and 
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FIG. 1: (a) Dispersion curves of a free atom for Ifc^l = 27r-4.5- 
W^m'^ and Q = 200s"\ (b) Density distributions |(/>o,iP for 
kd^2nx 4.5 ■ 10^m-\ fl = llOOs-S = u}± = 2-k ■ lOOs"^ 
and three different boson numbers: — 1390 (solid), 1420 
(dashed), 1460 (dot-dashed). 



radial (a;||) trap frequencies. As it results, this interac- 
tion clusters the majority of atoms within a single well, 
thus contrasting the quantum tunneling across the bar- 
rier. Due to the geometry of the problem, and taking 
into account that the scattering lengths aij are practi- 
cally equal, there is another ground state which is ob- 
tained from that of Fig. ^ by inverting the horizontal 
axis and interchanging 0o and 0i . Thus we have two sta- 
ble stationary states with equal energy. The numerical 
evidence of Fig. |^ is also supported by a synthetic vari- 
ational argument, already exploited in Ref. ^ for Li, and 
based upon a suitable two mode approximation. 

The matter field fluctuations allow the passage from 
one to the other state. If the decoherence is negligi- 
ble, coherent oscillations between such states can be ob- 
served, demonstrating MQC. To evaluate the oscillation 
frequency we quantize the two mode system, as done in 
Ref. ||. First of all, we write the classical Hamiltonian 
corresponding to the equations of motions (^, taking 
500 = 511 =901=9 



n 



ih^ka 
2m 



V'o*i?i0o+0rff/^i + |(i0op + i0ii 



(0*v0o - ri^A) - mri^o + roA) 



(fx (5) 



(from now on we omit the tilde on the 0's, even though 
we are always in the gauge Eq. (|^)). We then introduce 
the spinorial ground states 



V'Os(^) 
V'ls(^) 



; V's,2(^) 



(6) 



where -003 and ipig are the ground state wave-functions 
associated with the two internal states. 0g,2 is obtained 
from ijjg^i by interchanging the spinorial components and 
inverting the axes. 
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FIG. 2: First and second excited energy levels of the rubidium 
BEG versus the excess of atoms above the breakup value Ni 
for ti;|| — 27r • 100s~^ and some values of lo±. In the inset we 
plot Nd for the same values of the trap frequencies. 



FIG. 3: Three-body and two-body (inset) loss rates. The 
latter one refers to atoms trapped in the F — l,m — —1 level 
and represents an upper limit M3i. 



It is convenient to use the basis vectors V'a = fpg,i+'4'g,2 
and ipb = '03,1 ~ "03, 2- We write the quantized spino- 
rial field of the bicondensate as ipix) = aipaix) + bipbix), 
and substitute this expression in the operator version of 
Eq. (H). Thus, we obtain a reduced Hamiltonian Tired 
of the same form found in the attractive case of Ref. ||. 
By Ti-red, we evaluate the difference between the low- 
est eigenvalues, which provides the tunneling rate. In 
Fig. y we plot the first and second energy level as a 
function of the number of atoms and for some values 
of uj±. The energy of the ground state is set to zero, 
so the energy of the first level gives the tunneling fre- 
quency. The entanglement between the condensate and 
the lost atoms induces a decoherence of the superposi- 
tion. Using the approach of Ref. ||, we find that the 
coherence is given by C = e~*^*^, where M is the num- 
ber of lost atoms and e = 1 — 2/ cfix^ljQ{x)ipi{—x)/N. 
The quantity Nd = 1/e is the number of atoms which 
must be lost in order to reduce the coherence by 1/e. 
The inset of Fig. || shows how Nd scales with N. The 
relevant loss processes are two-body inelastic and three- 
body coUisional decays. In Fig. ^ we report the average 
three-body and two-body (inset) loss rates, the latter one 
refers to atoms in the m = — 1 level. We have used the 
upper limit of 1.6- 10~^^cm^ /s for the two-body loss rate 
coefficient and 5.8 • 10~^°cm^/s for the three-body pro- 
cesses, both of them measured in Ref. |l^ for the trapped 
Zeeman level F = l,m = —1. Two-body decay can 
occur also by means of collisions between atoms in dif- 
ferent Zeeman levels. We suppose that the correspond- 
ing loss rate is of the same order of magnitude as the 
measured value. From Figs. ^,|| we find that the deco- 
herence effects are negligible during a MQC oscillation 
period. For high uj± the threshold Ni and the loss rates 
are lower. For ojj_ — 2tt ■ 200 400s~^ the overall loss 
rate is much smaller than the corresponding tunneling 
frequencies. Therefore we can observe many oscillations 
before a single atom is lost. For u!± = 27r- 100s~^ the loss 
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FIG. 4: Distribution P{m) for some parameters. The solid, 
dashed and dashed-dot lines correspond to P{m) at the initial 
time, at a quarter of period and half a period, respectively. 
uj± is (a-b) 27r-100s"S(c-d) 27r-200s"^ and (e-f) 27r-400s"\ 

rate is > 5s^^. With a tunneling frequency of 8Hz also 
in this case we can observe an oscillation before a single 
atom is lost. If the loss of atoms does not transfer energy 
to the trapped atoms, the escape of a few atoms does 
not reduce the superposition coherence (inset of F igJ2 ), 
but modifies slightly the tunneling rate. In Ref. p^it 
is shown that inelastic coUisional processes induce local 
variations of the mean-field interparticle interaction and 
are accompanied by the creation/annihilation of elemen- 
tary excitations. This phenomenon depends on the den- 
sity and is completely negligible in our case. Notice that 
Eqs. ^ are invariant if N is varied by a factor a and the 
lengths and the energies are multiplied by a and a~^, 
respectively. So for uj = uj± = 27r • 70s"i the threshold 
is Ni = 1980. By the two mode model we find that the 
tunneling frequencies are reduced by nearly a factor 0.5. 
However the two-body and three-body decay rates are 
reduced by the factors 0.25 and 0.125, respectively. 

As we have shown in Ref. |2| the probability distribu- 
tion P{m) associated with the observable M = a^b + aw 
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displays two peaks for the ground state, that are the ahve 
and dead state of the Schrodinger cat. This observable 
is associated with the population measurement in one 
of the two Zeeman levels. To observe the coherence we 
have to put the system in one of the two states, that 
corresponds to take a superposition of the ground and 
first excited state of the condensate. This can be ob- 
tained measuring with a nondemolition technique (e.g. 
a phase contrast technique) the population in a Zeeman 
level when only the ground state is populated. This ob- 
servation collapses the P{m) distribution to one peak, 
as discussed in Ref. If the energy transfer is not too 
large we expect that only the first two energy states are 
populated. With this initial preparation the system be- 
gins to oscillate at the frequencies of Fig. |^ between the 
two Schrodinger cat states, as reported in Fig. ^. At the 
initial time only one peak is present (solid line). At a 
quarter of period corresponding to the frequency sepa- 
ration between ground and first excited state, the P(m) 
displays two peaks (dashed line). At half a period the 
only peak is that absent at the initial time, thus there is 
a coherent oscillation between the two states. Detecting 
such an oscillation would provide evidence of a SC at an 
intermediate time when both peaks are present. 

Notice that if we had N loosely coupled or independent 
atoms {Ng sufficiently low or even zero) the superposition 
of ground and first excited state would have a single peak, 
oscillating as a coherent state inside a harmonic poten- 
tial. This would by no means be a SC. On the contrary, 
we have shown that, for Ng sufficiently high, we have a 
two peak distribution with the two partial barycenters at 
nearly constant positions. During the evolution, the two 
peak amplitudes oscillate, that is, the probability to find 



the system in either state oscillates. 

In conclusion, we have shown that a double well dis- 
persion curve can be obtained by a suitable Raman cou- 
pling. In this situation a symmetry breaking in momen- 
tum space is demonstrated solving two coupled Gross- 
Pitaevskii equations. The condensate can oscillate be- 
tween the two emerging asymmetrical steady solutions 
(SC states) by means of the field quantum fluctuations 
(MQC). We have found that it is possible to obtain 
an oscillation frequency between the SC states around 
50 lOOs^^. In order to neglect the thermal activation, 
the second excited level energy {E2) has to be higher 
than the thermal energy. From Fig. ^we can see that E2 
ranges between 100 300s~^, that correspond to a tem- 
perature of 0.7^2.3nK . If the cooling is performed below 
the threshold, when the symmetry breaking does not oc- 
cur, the required temperature can be ^ 5nK. However 
it may not be necessary to cool at very low temperatures 
the whole condensate, but just the involved degrees of 
freedom, provided that this latter ones are weakly cou- 
pled with the other modes, which act as a thermal bath. 
We remark that for low densities the evaporative cooling 
allows to reach much lower temperatures, because of a 
smaller three-body decay rate. To be sure that no exci- 
tation is present one could tailor the trapping potential 
in such a way that only the first two levels are bound. 
In this work we have chosen the parameters for which 
the tunneling frequency is much larger than the deco- 
herence rate, however the symmetry breaking and the 
super-Poissonian atom fluctuations below threshold can 
be observed with a much higher number of atoms, thus 
these phenomena are observable with the present tech- 
nology. 
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